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Abstract 



The existing estimation of the Electric Dipole Moment (EDM) of the neutron in the 
SUSY SU(5) GUT model was quite smaller than one in the MSSM, despite of its 
increased degree of freedom. This paradoxical result can be resolved by appropriately 
estimating the CP-violation in the Higgs sector. As a result, the neutron EDM in 
this model is estimated lager than the existing one as expected. 

1. Introduction 

The Electric Dipole Moment (EDM) of an elementary particle is generated by the violation 
of CP0. In the standard model0,@, the neutron EDM (d n ) is theoretically evaluated as: 




(1) 



On the other hand, the experimental upper limit on the neutron EDM is: 




(2) 
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at present!!. However, this experimental limit will be improved in the near future. In several 
extensions of the standard model, it is well-known that the neutron EDM is expected or 
evaluated to be larger than one in the original standard model. Thus, once the neutron EDM 
is observed larger than the prediction Eq.(|l|) in the standard model, the reality of the other 
models is increased. In this paper, we estimate the neutron EDM in the SUSY SU(5) GUT 
model as one of the extensions of the standard model. Prior to the authors, the neutron 
EDM in this model is already evaluated by A. Romanino and A. StrumiaB. However, in 
their existing estimation^, there is not so much difference of the neutron EDMs between in 
the SUSY SU(5) GUT model and in the standard model except the large tan/3 settings, 
because the CP violating phase only in the quark sector of the mass matrix is taken into 
account, and such phase in the Higgs sector is essentially neglected. Generally speaking, 
according to the Kobayashi-Maskawa theoremS, a mass matrix consisted of 3 x 3 or more 
components can include the complex phase to violate CP. Thus, we introduce the complex 
phase into the Higgs sector of the mass matrix, because the Higgs/Higgsino superfields have 
five components in this model, and the existence of this new phase is one of the features of 
the SUSY SU(5) GUT model in contrast to the standard model or MSSM. When trying to 
estimate the neutron EDM in the SUSY SU(5) GUT model, this new mechanism to violate 
CP should be taken into consideration. In the following sections, all the estimations are 
carried out with including this new mechanism to violate CP symmetry. 

2. Estimation of Neutron EDM 

In this section, we present a sequence of the derivation from the general theory of SUSY 
SU(5) GUT to the estimation of the neutron EDM. 

A. Superpotential 

In the SUSY SU(5) GUT model, the superpotential is given as followJP: 
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Vs = Xi (jSP&Z + imE^Ej ) (3) 
+A 2 # a (Eg + 3m'^) if 6 

I / _ ttu yvw yah 

"T J jk^-uvwab^ 2 -A-j -<*-k 

ab\ 



+g jk H a X^Y h 



+h.c. , 

where, a and 6 are SU(5)'s, j and k are flavors' indices, respectively. On the Higgs super- 
multiplets H and H, H a is a 5-representation and ff a is a 5-representation. On the matter 
supermultiplets X and y, Xj b are the representation of 10 's, and Yj a are the 5's. Note that 

J£&b X^ a 

The complex phase parameter 9 to violate CP will be introduced into the second term 
of the superpotential Eq. (|3|) as the Higgs sector: 

3X 2 m'HH — ► 3X 2 m'He w H . (4) 
For more detail on the phase introduction, see the section 3. 

B. How to break SUSY SU(5) down 

It is well-known that most of the symmetries in this model, like SUSY, are broken in our 
real world; the sequence to break the symmetries down is as follows: 

1. The original SU(5) symmetry is breaking to SU(3)xSU(2)xU(l) with the energy scale 
getting lower by the Higgs mechanism of the heavy Higgs bosons in 24-representation. 

2. Moreover, the supersymmetry is broken spontaneously by added the soft breaking 
termJ as: 

mififa + (m A W 3 + Bono^hv + h.c.) + m 1/2 X j \ j , (5) 
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where, 0's are the scalar fields, W 3 stands for the trilinear scalar terms, hi and h 2 
are the two Higgs multiplets, and A's are the gauge fermions. m is the universal 
scalar mass, A is the universal trilinear coupling constant, B is the universal bilinear 
coupling constant, and mi/ 2 is the universal gaugino mass. On the meanings or the 
values of these introduced parameters are discussed in the section 3. 

3. Furthermore, SU(2) symmetry is spontaneously breaking to the electromagnetic gauge 
symmetry as usual: 



The observations of the neutron EDM are hold in this energy region. 
C. EDM formulae 

In order to estimate the EDM of a quark current, the shapes of one- loop graphs shown in 
Fig.l are regarded. The representations of each particles coincide with their physical states, 
where the explicit forms of the representations are shown in Appendix B. The diagrams 
including the virtual loop of the chargino, neutralino or the gluino can generate the CP-odd 
contributions to the current amplitude. Therefore, the contributions of these three kinds of 
the diagrams are taken into account. The graphs including the virtual loop of the heavy 
gauginos (i.e., X-ino or Y-ino) and the heavy Higgsinos are also possible to contribute CP- 
oddly, however, they are neglected. Roughly speaking, since the EDMs are proportional to 
the inverse of the virtual particle masses ; the EDM values of such graphs including heavy 
virtual particles are extremely small compared to the graphs of the charginos, neutralinos, 
or the gluino. 

For the purpose of calculating the loop correction to a current of the quark with spin 1/2, 
the current amplitude is decomposed according to the form factors as follows: 



SU(2) x U(l) -> U(l) 



em 



(6) 



{qfinit(p')\f(q)\qinitial(p)) 



(7) 



4 



Qfinitip') 



2m 



+ 75^ qu 



2m 



-\^-q»)isF A tf) + 



qinitial 



where, j^(q) is the current with four- momentum transfer q = p' — p. F 3 (q 2 ) is the CP-odd 
form factor. 



The EDM(d/) of a quark with the flavor / is given by 

df = ~^ F M 2 - 0) . (8) 

By making use of this definition, the neutron EDM(d n ) is written as: 
4 1 

d n = ~dd - ~d u , (9) 
where, dd means the d-quark EDM, as d u is the u-quark EDM. 

Moreover, the Chromoelectric Dipole Moment (CDM)0,0is also taken into account. The 
CDM is defined as the factor d CDM of the effective operator arise in the QCD Laglangian: 

Ccdm = - l -d CDM qa^T a qG^ a , (10) 

where, T a, s stand for the generators of SU(3). 

The CDM contributions to the neutron EDM0 is as follows: 

dn = \e [\dT M + \<f u DM ) ■ (11) 
For more details on the loop calculations of the EDM and CDM formulae, see Appendix C. 



3. Evaluation of parameters 

In order to evaluate the soft breaking mass parameters, on the correction of the renormal- 
ization group flow at one loop level is discussed in this section. The explicit forms of the 
Renormalization Group Equations(RGEs) are shown in Appendix A. In order to solve the 
RGEs, we assume the following conditions. 
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The gauge coupling constants are bent at Msusy — l(TeV), and unified at M S u{b) — 
10 16 (GeV)§. 

Yukawa coupling constants of the top and bottom quarks are given to coincide with 
their values at electroweak scale Mz independent of tan/3. 

The spontaneous SUSY breaking down scale M x is as: 

M x = Gp 1/A Mp /2 ~ 6.0 x 10 10 (GeV) , (12) 

where, Gp is the Fermi coupling constant and Mp is the Planck mass scale. Therefore, 
Mx is fixed to 6.0 x 10 10 (GeV). Note that this condition is ad-hoc for the pure SUSY 
SU(5) GUT, because this is a consequence of the super Higgs effect in hidden sector 
based on N=l supergravity grand unified theory!. Also this is mere a formal boundary 
condition, indeed the values of and Aq are almost independent of Mx- 

All the masses of the supersymmetric scalar particles are set equal to m at energy 
scale Mx- 

Bq = (Aq — l)m is assumed. This is another ad-hoc condition for the pure SUSY 
SU(5) GUT adopted from the minimal supergravity theories!. 

Once an assumed Bq at Mx is given, the parameters like mhi or m/^mass parameter 
of two Higgs doublet) at Mz scale can be obtained by solving the RGEs directly. On 
the other hand, to break the electro-weak symmetry down as usual, two conditions: 



= -(m 2 hl + m 2 h2 + 2/i 2 ) sin 2(3 , (13) 
, 2 _ m l , ™? hl - m\ 2 tan 2 (3 

^ ~ 2 + tan 2 /3-l ' 1 j 

are required. As the results of these two equations, B and fi at M z are derived, 
however, this value of B at Mz should be coincident with another estimation. To 
match both two B values at Mz, the initial Bq at Mx is tuned recursively and decided 
by numerical analysis. Thus, Bq ( and also Aq) is fixed. 
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vii) tan (3 is restricted as 2 < tan j3 < 40, to avoid emerging the Landau pole divergence of 
Yukawa coupling constant of the top quark. 

viii) The other parameters are cited from the last report of Particle Data Group0. 



As it is slight confusing because there exist many introduced parameters with mass dimen- 
sion, we give a summary of them in the Table 1. 

With these conditions, several additional restrictions on other parameters are obtained 
from the RGE analysis. The existence of m.1/2 in Eq.(j|) implies the universality of all 
the gaugino masses at the SUSY breaking energy scale Mx = 6.0 x 10 10 (GeV). The condi- 
tion Msu(5) > Mx renders this universality, because, the unification of the gauge coupling 
constants is broken at M su ^ in the first, however, the supersymmetry to make the gaugi- 
nos massless still remains until the energy Mx- As the results, the following three relations 
between mi/2 and each gaugino masses or gauge couplings are realized at one loop leveM 
mi(M z ) _ m x /2 



(2 = 1,2,3), (15) 



oii(M z ) OiiMx) 

where a^'s are the gauge coupling constants. Note that these relations give each (different) 
gaugino masses at Mz to solve the RGEs, however, they never mean the universality among 
the gauge couplings at Mz or Mx, i.e., a* 7^ atj(i 7^ j), in general. 

Furthermore, several constraints are derived among A2, mf, and the CP-violating phase 9' 
in the Higgs sector. The parameter [i can be regarded as the Higgsino mass, and is defined 
as: 

/i = A 2 (3m - 3mV e ) . (16) 

By our RGE analysis, is about l(TeV). As m implies the VEV of the Higgs bosons in 
the 24-representation, thus the order of m is 10 16 (GeV). The unified mass of the 5 and 5's 
SU(3) heavy Higgs bosons (m co i or ) is given as: 

m color = X 2 (2m + 3m'e w ) . (17) 



3 1 Tricolor I 


- 2 




3(|m coior 


l + H) 



Quite small contributions to m co \ OT from the soft breaking terms Eq.fl5|) are neglected. The 
order of the mass m co i or is estimated as 10 16 (GeV) by the proton decay analysis!!. Due to 
these two relations Eq. (|T6|) and Eq. (|T7|) , the phase 6 is restricted as: 

| sing] < sinfl max = - 5 l m ^^| „ 0(10-13) (18) 

and the ratio m'/m is restricted as: 

< ^ < 3\ m color \ + 2|/i| , , 

_ m ~~ 3(|m co , or | - 

with the quantitative relation between |m co / or | and 

\m co i or \ » . (20) 
Therefore, 

m'/m ~ 1 , (21) 

is valid in good approximation. This result coincides with the fact that m! is introduced as 
the fine-tuning parameter to keep the Higgsino mass in the usual region ~ l(TeV), i.e., the 
order of m! should be the same with m as 10 16 (GeV). On the other hand, the parameter A2 
is written as: 

A — \ m color\ (22) 

2 ~ \2m + 3m'e id \ ' 1 ' 

therefore, the order of the parameter A2 is evaluated as O(10 _1 ). When /i is rewritten as: 

[i = |/i| exp {16') , (23) 

the complex phase 6' of /i is given by: 

sine' = -r-rm'lsinfll . (24) 

H 

The phase 6' of /i as the CP-violating phase in the Higgs sector is given by this function 
of A2, m' and sin^. Obviously, the maximum sin# as O(10 -13 ) corresponds to the 
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maximum CP- violating phase 9'. The order of the maximum sin#' is derived from Eq.(p4|) 

as: 

| sin 6*| =^ sm9 max , sin#' ~ 1 , (25) 

by making use of two conditions Eq.([l8]) and Eq . (^2|) . In Fig. 2, the phase dependence of the 
neutron EDM is given as a function of sin#. In the following sections, 9 = 9 max is fixed. 

4. Numerical Analysis and Results 

The rest free parameters are m^, m , and tan/?. These parameters are expected to settle 
in such regions as usual: 

• m 2 : TeV order. 

• m : less than a few TeV . 

• tan/? : 2 ~ 40 (~ \m t \/\m b \ at M z ) 

In Fig. 3, the neutron EDMs are shown as functions of raxn : 1 < mj/2 < 10 (TeV) for 
tan/3=3, 5, and 10. mo =l(TeV) is fixed. The flat line is the experimental upper limit of 
EDM. Obviously, m^'s value is required to be larger than 4.5(TeV) by the experimental 
limit. 

In Fig. 4, on the m\/2 dependence of the neutron EDM is drawn for mo =1, 3, and 5 (TeV). 
tan/3 = 5 is fixed. The flat line means the experimental upper limit as in Fig. 3, min is 
necessary to be larger than 5 (TeV). The curves imply that the neutron EDM is not sensitive 
on m . 

Figure 5 summarizes the topographic plots of EDM's experimental upper limit as func- 
tions of mo and mi/2 for tan/3=3, 5, and 10. These plots show the two facts - Toil's value 
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should be larger than a few TeV, and tan (3 is required to be small. The EDM is more 
sensitive on m-1/2 rather than m . 

5. Conclusion and Summary 

In summary, the neutron EDM in the SUSY SU(5) GUT model is obtained by introducing 
the complex phase into the mass matrix of the Higgs sector. As the results, the neutron 
EDM of this mechanism is derived significantly larger than existing estimations by the CP 
violation in the quark sector. According to this results, the constraints on the m , ^1/2 
and tan f3 become more strict than existing evaluations, however, the allowed parameter 
region still remains. Additionally, we should take care of the phases of other parameters, 
like m^j's, because such phases can cancelEl the effect of the introduced phase into the 
Higgs sector. If such phases were possible to exist, the evaluation of the EDM in this paper 
should be reconsidered as the maximum, i.e., the strictest estimation. Anyway, once the 
neutron EDM would be observed larger than the prediction of the standard model Eq. ([[]), 
the other models beyond the standard model, including the SUSY SU(5) GUTs(predicting 
large EDM), will be more realistic. 
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Appendix A: Equations of Renormalization Group 

The RGEs at one loop level are as follow|l,i: 
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da a 

~dt 

dm a 

~df 

da t 

~dt 
dab 
~~dt 
d_M 
dt 
dAu^c 

dt 
dAr 





(Al) 


-2b a a a m a , 


(A2) 


2a t {-c K a >a a + 6a t + a b ) , 


(A3) 


2a b {-c\ >a a + a t + Qa b ) , 


(A4) 


-2c^a a m a /m + 12a t A i + 2a b A~ b , 


(A5) 


-2c^a a m a /m + Qa t A~ t , 


(A6) 


-2c^a a m a /mo + 2a t A~ t + 12a b A~ b , 


(A7) 


-2c^a a m a /m + 6a b A~ b , 


(A8) 


—2c\ 'a a m a + oa t m Af + ba b m A~ b , 


(A9) 


— 9r^rv m 2 4- RnivY 2 


( Ami 


-2c (3) a„m 2 + 6m£ 2 , 


(All) 


-2c (4) a„m 2 + -aiTr(Ym 2 ) , 


(A12) 


-2ci 5) a a m 2 - -aiTrfFm 2 ) , 

CL C* Q, p- y ? 




(A13) 


2 

— 2r^(y m 2 A n^r(Ym 2 \ 

Lx a lll 'a' r Lx l 1L \ 1 "* ) i 






2a t E; + 2a b i: b - 2c^'a a m z a + -aiTr(Fm^) , 




(A15) 


AatT 2 - 2cfa a m 2 a - ^a 1 Tr(Fm 2 ) , 


(A16) 


4« 6 S 2 - 2cfa a m 2 a + ^a 1 Tr(Fm 2 ) , 




(A17) 



dt 

dt 
dB_ 

~dt 

dm 2 hl 

dt 
dm 2 h2 

dt 

dm 2 _ j - r 

u,c,d,sL 

eft - 
eft 

dt 

dm b L 
dt 

dm iR 

dt 

dm lR 

dt 



where, the dimensionless energy scale variable t is defined as t = (l/47r) ha(Q/Msu(S)) derived 
from the bare energy scale variable Q. m a 's are the gaugino masses, and a a = g 2 /47r's are 
the gauge coupling constants. Ag and trig are squarks' trilinear coupling constants and their 
masses, respectively. The suffix L or R indicates each partner's chirality. The suffix a is 
always summed up from 1 to 3. The optional coefficients and variables are defined as: 
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, -41/10, 19/6, 7) , M Z <Q< Msusy , 
b a = { , (A18) 

^ ;-33/5, -1, 3) , Msusy <Q< M su{5) 

c« = (13/15,3,16/3), (A19) 

ci 2 ) = (7/15,3,16/3), (A20) 

cf = (3/5,3,0), (A21) 

#> = (1/15, 3, 16/3), (A22) 

= (16/15,0,16/3) , (A23) 

cf = (4/15, 0,16/3), (A24) 

E 2 = (Afml + m 2 2 + m\ L + m 2 iR ) , (A25) 

S 2 = (A\ml + m 2 hl + ml + m\ R ) . (A26) 

All the Yukawa couplings except of the top and bottom quarks are neglected, because their 
values are too small to be taken into account. The exceptional two (relatively large) Yukawa 
couplings a t f) = gtb/^ 71 are defined as usual: 
92 m t 1 

9t = -7= > ( A27 ) 

V 2 m\v sin p 

^ = -75 a • A28 

V2 cosp 

Furthermore, keeping the traceless conditions of the SU(5) gauge: 

Tr(Ym 2 ) = m 2 Tr(Y) = , (A29) 
where, the 5x5 diagonal matrix Y is defined as: 

Y = diag{-2/3, -2/3, -2/3, 1, 1) , (A30) 
are required to avoid the gravitational mixed anomalyEl. 

Since the initial or boundary conditions are necessary to solve the RGEs, they are given at 
the electro-weak scale Mz — 0.9 x 10 2 (GeV) for the gauge and Yukawa couplings, and at 
Mx — 6.0 x 10 10 (GeV) scale for the other variables. Note that we set all the initial values of 
the soft breaking mass parameters equal to mo, and an ad- hoc constraint Bq = (Aq — l)mo 
from the minimal supergravity theory is applied. 
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Appendix B: Representations of particle states 

As one of the general features of the quantum theory, the mixture among the particles with 
the same quantum numbers usually occurs, then, some non-diagonal elements arise in the 
mass matrix of bare particles. In order to coincide their representations with the physical 
particles, the mass matrix should be diagonalized and realized as possible as they can be. 



The gauginos and Higgsinos are mixed as the charginos and the neutralinos. As the results 
of the mixture, the mass matrix of the bare charginos is : 

m 2 \f2mw sin (3 
\[2m w cos (3 X 2 b 



Mr 



( 

\ 



(Bl) 



and of the bare neutralinos is 
/ 



M 



N 



V 



/ 



(B2) 



mi — m z sin 9 W cos (3 m z sin 9 W sin (3 

m 2 mz cos 6w cos (3 — mz cos 9w sin (3 

— mz sin 6 W cos (3 m z cos 9w cos j3 — X 2 b 

mzsinOwsmP — mz cos 9w sin f3 — X 2 b 
where, m\ and m 2 stand for the gaugino masses, and mz are W and Z boson masses, 
respectively, b = 3m — 3mV e ; m is the Vacuum Expectation Value(VEV) of heavy Higgs 
bosons in the 24- representation, tan (3 is defined as tan/? = v 2 /vi; v\ and v 2 are the VEVs 
of the Higgs scalar fields 5 and 5, respectively. These mass matrices are diagonalized by the 
biunitary transformation for the charginos! as: 



U*M C V = diag(\m Cl \, \m C2 \) 



(B3) 



where, the suffices Ciii = 1,2) represent the physical chargino states. Moreover, the neu- 
tralinos are diagonalized by the complex orthogonal transformation as: 



N t MnN = diagQm^l, |mjv 2 |, |mw 3 |, l 77 ^!) 



(B4) 



where, the suffices iVj(i = 1 ~ 4) imply the physical neutralino states. Note that the repre- 
sentations of these diagonalized physical states as the charginos and neutralinos are used all 
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through this paper, instead of the bare fields of the Higgsino and gaugino. It is not necessary 
for the gluino to diagonalize the mass matrix in contrast to the charginos or neutralinos. 
Therefore, the gluino mass m- g is simply equivalent to in Eq.(|l|). We introduce a symbol 
g as the gluino. 



The bare mass matrix (6 x 6) of the u-type squarks is given as: 

m\ L + D L 1 + mjj m u (Atm + X 2 b(v l /v 2 )) ^ 
mu{Afn Q + X 2 b*(v 1 /v 2 )) m\ R + D R 1 + m 2 v 

and of the d-type squarks is: 



M, 



u 



(B5) 



M 



D 



m 2 - L +D L l + m 2 D m D (A*-m + X 2 b(v 2 /v x )) 
m D (Afn + A 2 &*(i>2A>i)) m\ R + D R 1 + m 2 D J 



(B6) 



where, 



D L = m 2 z (T 3UiD - Qu,d sin 2 9 W ) cos 2/5 
Dr = m 2 z Qu,D sin 2 9w cos 2(3 . 



(B7) 
(B8) 



Qu,d stands for the electric charge of the u, d-type quark, respectively. T^u = +1/2, and 
T 3 d = —1/2. mjj,D denotes the u, d-type ordinal quark mass matrix(3 x 3), respectively. 
rriqi and rriq R are the matrices of the soft breaking mass parameters of the squarks. Aq 
is the trilinear coupling constants matrix. Both of the squark mass matrices(6 x 6) are 
diagonalized as: 







l' m «21 1 1 1 


l,K~J, 


l m d 13 l> 


1 £121 1 ' 1 



\ m U23 |) 



(B9) 
(BIO) 



where, the suffices u, (d)kj(k = 1,2) mean the physical states of the u,(d)-type squarks of 
the j-th generation, respectively(the sub-index k stems from two chiralities of the ordinal 
quarks as the superpartners of the squarks, however, this has no physical meaning). These 
definitions of the squark states are used all through this paper, like of the charginos and 
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neutralinos. Note that we neglect all the non-diagonal elements of Aq, rriqL,R and mu,D and 
their additional phase in the squark mass matrices to violate CP, however, this additional 
effect has already estimated by others§,i; the EDM generated by this additional phase is 
quantitatively quite smaller than one by the phase of the Higgs sector in our discussions. 



Appendix C: Neutron EDM Formulae 



The neutron EDM and CDM are given by following equations^!,!. All the notations of the 
transformation matrices and their elements, like V* 2 in Eq . (|C5p , are the same with Appendix 



B. The chargino contribution is: 

2 2 

,C I a em \ -> \ T / p O \ m Ci 



d)/e 



47rsin 6 W 



2 

m- 



f'k 



Q f ,B 



m Ci \ , in n \ a i m Ci 



I 2 



+ (Q f -Qf,)A 



jC-CDM 



f'k, 
2 2 



in 



f'k, 



47rsin 6 W j^i^J 



m 



a 



fikJ 2 I 7r) 2 I 

f'k \ f'k, 



(CI) 

(C2) 



where, f = u,d for /' = d, u, respectively. a em is the fine-structure constant, g% = \/Attos, 
and 9w is the Winberg angle. Each definitions of two introduced functions A(r) and B(r) 



are: 
A(r 

B(r) 



1 



2(1 -r) 2 
1 



21nr\ 
3-r + - ) , 



1 + r + 



1 — r J 
2r lnr 



2(1 -r) 2 V ' ' 
respectively. The T's are defined as: 

^uik = K uVi2Ddlk(U* 1 D^ lk — ^dU* 2 D* d2k ) , 

F%k = K d U* 2 D ulk (V*D* ulk - K u V*D* u2k ) , 
and 

m,. 



\/2mw sin (3 ' 
\/2m w cos j3 



(C3) 
(C4) 

(C5) 
(C6) 

(C7) 
(C8) 
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"I 














9 \ 


f m- Ni 


i m 2 r. / 







The neutralino contribution is: 

Tf ifc = [- v / 2{tan^(g / - r 3/ )JV u + T 3f N 2i }D* flk + K f N bt D} 2k ] 

x (v^2 tan 6 w Q f N H D f2k - K f N bi D m ) , (Cll) 

where, the suffix b = 3(4) for T 3 f = — l/2(+l/2), respectively. 
The gluino contribution is: 

4/--^|>(ry^^), (Ci2) 

4- CTi '^gi m (ry^(i|) (c 13) 

where, 

T% = D f2k D* flk , (C14) 

w-ifrM 1 "^-^)- (C15) 

Note that all the generations of the squark are summed up in the virtual loops. At last, in 
order to compare with the EDM of the real neutron, the total sum of above EDM and CDM 
contributions is composed subject to the formulae Eq.(|D and Eq.flTTD. 
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TABLES 



Table 1. Introduced parameters with mass dimension 



Symbols 


Energy scales 


Descriptions 


M P 


10 19 (GeV) 


The Planck mass scale. 


M SU(5) 


10 16 (GeV) 


The unification scale among the gauge coupling constants. 


M x 


6.0 x 10 10 (GeV) 


The spontaneous supersymmetry breaking down scale. 


MsiJSY 


l(TeV) 


The energy scale emerging the effects of the gauginos. 


M z 


91(GeV) 


The electro-weak scale 
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